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Abstract 

We address numerical aspects of local quark-hadron duality using the exam- 
ple of the exactly solvable 't Hooft model, two-dimensional QCD with N^. — > oo. 
The primary focus of these studies is total semileptonic decay widths relevant 
for extracting \Vcb\ and \ Vub\- We compare the exact channel-by-channel sum 
of exclusive modes to the corresponding rates obtained in the standard I/wq 
expansion arising from the Operator Product Expansion. An impressive agree- 
ment sets in unexpectedly early, immediately after the threshold for the first 
hadronic excitation in the final state. Yet even at higher energy release it 
is possible to discern the seeds of duality-violating oscillations. We find the 
"Small Velocity" sum rules to be exceptionally well saturated already by the 
first excited state. We also obtain a convincing degree of duality in the differ- 
ential distributions and in an analogue of R^+e- (s)- Finally, we discuss possible 
lessons for semileptonic decays of actual heavy quarks in QCD. 
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1 Introduction 

Questions of how to compare hadronic observables to the apparent underlying fun- 
damental theory of QCD lie at the heart of understanding the nature of strong in- 
teractions. Thirty years after its inception, QCD in D = 4 spacetime dimensions still 
stubbornly refuses to admit a global solution. The asymptotic freedom property of 
the theory permits the perturbative calculation of (Euclidean) Green functions in- 
volving large values of momentum transfer or energy release in terms of quarks and 
gluons, the fundamental objects of QCD. But at lower scales one enters the nonper- 
turbative regime, which not only invalidates (or at least complicates) the standard 
perturbative methods of field theory developed in QED, but also leads to a dramatic 
change in the physical spectrum of the theory. Instead of quarks and gluons, only 
colorless hadrons are produced as asymptotic states in any process, even at arbitrarily 
large energy. 

Many nontrivial theoretical techniques respecting QCD first principles have been 
developed to study nonperturbative features of the theory. Yet despite numerous 
advances, no one has been able to compute the masses, wavef unctions, or transition 
amplitudes of hadrons in terms of quark masses and couplings directly from the QCD 
Lagrangian. Moreover, many existing theoretical tools are expressed through various 
expansions in certain small parameters; the actual range of each parameter where the 
expansions are applicable is often not well known. In such a situation, it is clearly 
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advantageous to build a soluble toy field theory that incorporates as many features 
of the QCD Lagrangian as possible. 

Such a theory does indeed exist, the famous 't Hooft model which is defined 
by the Yang-Mills Lagrangian in D = 2 spacetime dimensions in the limit of a large 
number of colors Nc- As was shown in the original paper, the quark-antiquark sector 
of the theory admits an infinite tower of confined, color-singlet solutions that can 
be obtained, in principle, to an arbitrary degree of numerical accuracy. The rea- 
son for this solubility lies precisely in the defining features of the model. Large Nc 
eliminates all Feynman diagrams with internal qq loops and nonplanar gluons. On 
the other hand, D = 2 allows gluon self- couplings to be eliminated by gauging away 
one component of the gauge potential A^. Since only two components are initially 
present, the commutator term [A^, A^] in the covariant derivative, which gives gluon 
self-coupling, vanishes identically in such gauges. Then the only remaining Feynman 
diagrams to be summed for the quark-antiquark Green function are "rainbow" and 
"ladder" diagrams, whose Schwinger-Dyson equations can be solved, giving rise to 
an integral expression called the 't Hooft equation (discussed in Sec. 

The 't Hooft model provides an excellent laboratory for testing various approaches 
to strong interaction physics. After all, the 't Hooft equation provides a means to 
compute hadronic masses, wavefunctions, and transition amplitudes in terms of the 
underlying partonic degrees of freedom. 

In this work we are specifically interested in questions of local quark-hadron du- 
ality in the inclusive decays of heavy quarks. The notion of duality in general terms 
was first introduced in the early days of QCD in Ref. but not pursued for quite 
some time. A more detailed consideration was given a few years ago by Shifman 
and later reiterated in a number of papers (see, e.g., Refs. @, |^), with applications 
relevant to Minkowskian observables amenable to study via an operator product ex- 
pansion (OPE). This allows the formulation of the concept of local duality in a more 
quantitative way, including nontrivial nonperturbative effects; we refer the reader to 
these recent publications for the theoretical aspects. Here the question of duality 
is studied concretely by comparing the weak decay width of a meson containing a 
heavy quark computed in two ways. In terms of partonic degrees of freedom, one 
has an OPE depending upon the free quark diagram (with perturbative corrections) 
and a number of nonperturbative matrix elements suppressed by powers of the heavy 
quark mass. In terms of the hadronic degrees of freedom, one simply computes the 
weak decay amplitude for each allowed exclusive channel, and adds them up one by 
one. This comparison is especially instructive since one may consider the behavior of 
solution as the mass uiq of the heavy decaying quark is varied. 

Such a problem was first considered in Ref. 0, where the main elements in nu- 
merical computations of exclusive decay rates were annunciated. The hadronic result 
was compared to the Born-level free-partonic diagram as a function of niq. In terms 
of the OPE, the latter is the tree-level piece of the Wilson coefficient corresponding to 
the unit operator. The numerical agreement was seen to be remarkable, in that the 
onset of the asymptotic agreement was clearly visible already for relatively small val- 
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ues of niQ. The intrinsically limited numerical accuracy for sufficiently heavy quarks, 
however, prohibited drawing a definite conclusion about the size of nonperturbative 
corrections for asymptotically large mq. Additional numerical studies 0] consid- 
ered similar questions for weak decay topologies other than the simple spectator tree 
diagram, in particular weak annihilation (WA). 

The validity of the OPE was addressed analytically in Refs. 0, 0], which con- 
sidered on one hand the nature of the OPE for heavy quark decays, and on the other 
an explicit l/mg expansion of the decay amplitudes, which allows an analytical sum- 
mation of the individual decay rates in the asymptotic regime. The agreement of the 
two approaches through relative order l/mg was obtained by means of a number of 
sum rules derived directly from the 't Hooft equation, the archetype of which first 
appeared in Ref. |]TU| . 

While adequate to illustrate the theoretical validity of the OPE for the inclusive 
decay widths of heavy fiavors, the analytic methods per se cannot help in answer- 
ing the practical question relevant to phenomenology of beauty and charm quarks: 
Namely, how accurately do the OPE-improved parton computations describe the true 
weak decay width of a heavy fiavor meson with finite mass, only a few times larger 
than the typical strong interaction scale? A purely analytic expansion can hardly 
be used for this purpose, since it is a priori unknown how small an expansion pa- 
rameter must be for the expansion to start yielding a reasonable approximation, not 
to mention achieving the necessary precision. To obtain insights into the size of 
deviations between the actual decay widths and the expressions obtained from the 
OPE for quarks in the intermediate mass range, one must employ real numerical 
computations. 

In this paper we focus on semileptonic decays of heavy quarks. In the contexts of 
both real QCD and the 't Hooft model, they are technically simpler than nonleptonic 
decays. Moreover, the magnitude of local duality violation is phenomenologically 
most important in semileptonic decays when one extracts \Vcb\ and \Vub\- We use the 
techniques developed in Ref. to evaluate the required decay rates, and confront the 
total decay width with the expansion in terms of a power series in l/mg of Ref. 
Moreover, by making use of a number of relations derived in the large-mQ limit of 



the 't Hooft equation []TT|, members of the set of nonperturbative matrix elements 
involved can be related to each other, providing an economical description of the 
nonperturbative physics. These are the tools that allow us to study the onset of 
quark-hadron duality. 

As explained in Appendix |^, we use a scheme based on the modified Multhopp 
method, by which the 't Hooft equation is converted into an infinite-dimension eigen- 
vector system that for practical reasons must be truncated at some number N of 
eigenvector modes. The asymptotic convergence of this approach has not been rigor- 
ously studied, although it apparently must yield unlimited accuracy when the number 
of the Multhopp modes goes to infinity. Yet the rate of convergence at large 
is not well known. Additionally, large quark masses turn out to require one to use a 
larger A^ for sufficient numerical accuracy, as discussed in Sec. 0. It therefore seems 
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mandatory to make an independent cross-check of the numerical accuracy. We inves- 
tigate this problem by comparing the numerical values of a number of static properties 
of heavy mesons at different values of mq, with the results of their l/mq expansions 
obtained analytically from the 't Hooft equation; this is the topic of Sec. |[ We find 
that our solutions have sufficient numerical accuracy for masses rriQ corresponding to 
physical values (in the sense explained in Sec. |^) as large as 20 GeV. 

The duality of the inclusive widths of heavy-flavor hadrons to the parton-level 
widths, including the power corrections from the OPE, emerges through a set of sum 
rules that equate sums of weighted transition probabilities to possible final states 
and expectation values of the local heavy quark operators. Since our main interest 
lies in 6 — i> c transitions, which carry in practice a limited energy release, the most 
relevant are the so-called small velocity (SV) sum rules, which we study here in the 
heavy quark limit. The behavior of these sum rules not only shapes the semileptonic 
b ^ c decays in actual QCD, but is also important for the determination of the basic 
parameters of the heavy quark expansion. 

An additional advantage of the heavy quark limit for our investigation is that we 
are able to compute the SV amplitudes semi-analytically, using the exact relations 



Tl| derived from the 't Hooft equations and relying for input only on a few static 
parameters, which can be computed with a high precision. A discussion of these 
relations appears in Sec. We find that the SV sum rules in the 't Hooft model are 
saturated to an unexpectedly high degree by the first excitation above the ground 
state (which we henceforth call the "P-wave" excitation, despite the fact that in D = 2 
only radial excitations occur). Its contributions to even the Darwin {p%) and kinetic 
(/i^) expectation values constitute over 90% and 96% of the totals, respectively, while 
it saturates the "optical" sum rule for Mg — mj, to a 1.5% accuracy. This appears 
to be an intriguing dynamical feature of the model. A similar high-saturation effect 



has been observed in a quark fiux-tube model [|T^], for the contribution from the 
"valence" quarkonium states. 

We study the size of violations of local duality in the semileptonic decays b ciu 
assuming vectorlike weak currents and massless leptons. These assumptions are im- 
portant for comparison with QCD far beyond the obvious parallel of closely resem- 
bling the actual world: The strength of the resonance-related duality violation cru- 
cially depends on the threshold behavior in the decay probabilities, which is com- 
pletely different in two and four dimensions. The two-body phase space, while oc \p\ 
in D = 4, is oc that is, infinite at threshold, in D = 2. On the other hand, the 

situation is special for massless leptons: Their invariant mass is always zero if they 
are produced by a vectorlike source, and the weak vertex is then proportional to the 
momentum. As a result, in this case the threshold behavior of the decay rate becomes 
oc \p\ much in the same way as in real QCD. This is a crucial detail if one tries to 
draw practical lessons from the 't Hooft model. The need for a vectorlike coupling 
in D = 2 is even more stark for the parton-level calculation. There one finds that the 
integrated three-body phase space actually diverges for massless leptons, and only 
the behavior of the weak decay amplitude renders the width finite. We provide more 
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arguments in favor of such a choice in Sec. ^ which is dedicated to the inclusive decay 
widths. 

In Sec. ^ we briefly illustrate how well the duality works for the vacuum correlator 
of light quarks in the timelike domain. In the context of the heavy quark expansion 
this is relevant for the nonleptonic decay widths, including spectator-dependent effects 
hke WA. 

Section |^ summarizes our investigation and discusses the conclusions that can be 
drawn for actual QCD. 

Appendices describe the computational technique employed and contain a number 
of relations for the heavy quark limit of the 't Hooft equation employed in these 
numerical studies. 



2 The 't Hooft Equation and Its Solutions 

We first review some well-known properties of the 't Hooft model both as a reminder 
and to establish notation. Confinement is manifest in 1+1 spacetime dimensions with 
large Nc, and the quark(mi)-antiquark(m2) two-particle irreducible Green function, 
i.e., the meson wavefunction ip{x), is given by the 't Hooft equation: 



m 



M„^ ^^{x) = { + ) </.„(x) - (3' I dy ^n{y) Pt:— , (i) 



X 1—x J Jo [y—x 



where x is the momentum fraction in light-cone coordinates carried by the quark, 
and 

f3'^fiN,-im. (2) 

Since (3 is finite in the large-A^c limit, it provides a natural unit of mass. Thus, all 
masses in this paper are understood as multiples of /?. Indeed, as pointed out in 
Ref. P, P fills the role in 1+1 dimensions of served by Aqcd in 3+1. We discuss the 
estimation of /3 particular number in Sec. 0. 

The singularity of the QCD Coulomb interaction in Eq. is regularized using a 
principal value prescription, indicated by P in Eq. (|l]). 

Solutions n = 0, 1, ... of the 't Hooft equation alternate in parity, with the lowest 
being a pseudoscalar. The general analytic solution in closed form is not known. 
As the eigenvalue index n increases, the eigenvalues asymptotically approach 
(3'^[-iT^n + 0{\nn)]. 

The static limit mi = mg — >^ oo is most easily studied H, [l^, |ll| by employing the 
"nonrelativistic" variables M„ = mQ + e„, t = {l—x)mQ and \I'„(t) = -^^(^„(^1 — 
in terms of which Eq. (|1]) assumes the form 

+ f d,|^. (3) 
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We solve the finite-mass 't Hooft equation using a numerical method called the 
Multhopp technique |T^, a venerable system for solving integral equations with sin- 
gular kernels. It was first applied to the 't Hooft equation in Ref. |M. The idea 
is to expand the wavefunction in a series of modes, not unlike Fourier analysis, and 
then turn the equations for the mode coefficients into an equivalent infinite-dimension 
eigenvector problem. In practice, one then truncates at some point where the higher 
modes are deemed to have little effect upon the wavefunction solutions, which is of 
course strongly dependent on the highest value of n used. The detailed formulas for 
applying the standard Multhopp technique to mesons with unequal quark masses in 
the 't Hooft model appear in Appendix A of Ref. 0. 

Intrinsic to the original Multhopp technique is the evaluation of the wavefunction 
at a discrete set of points called "Multhopp angles," which in the current problem 
are equivalent to 



1 + cos 



fcvr 



N + 1 



k 



N 



(4) 



where is the number of modes retained in the numerical solution. The mode 
coefficients are then obtained by the use of a discrete inversion formula [(A7) in 0]. 
However, the Multhopp solutions can be seen to vanish as ^/x and y/l—x at the 
endpoints x = and x = 1, respectively [see (AlO)-(All) in 0], while the exact 
solutions are known to vanish as x'^^ and [l — x)"''^, respectively, where 



,2 + 7r7iCot7r7i 



(5) 



leading to a type of Gibbs phenomenon in the Multhopp solutions. Since the Mul- 
thopp angles cease to sample the wavefunction at some finite distance from the end- 
points, it may be expected that the wavefunctions thus obtained are numerically 



inaccurate there. This shortcoming led Brower, Spence, and Weis [0 to improve the 
Multhopp technique by eliminating the Multhopp angles and using instead a contin- 
uous inversion formula. The algebraic details are presented in Appendix ^ and it is 
this improved numerical technique that is used in obtaining our results. 



3 Heavy Quark Expansion and Cross Check of the 
Algorithm 

Let us first establish a bit of notation. The mass of a heavy quark of flavor Q is 
labeled as mi rriQ] in the weak transitions considered in subsequent sections, the 
final-state quark q is assigned the mass rriq. The spectator antiquark mass 1712 is 
labeled by m, or rrisp if there is any chance of confusion. 

As explained in the previous section and Appendix 0, we use the modified Mul- 
thopp technique to find numerical solutions of the 't Hooft eigenstate problem. Since 
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the heavy meson wavefunctions are peaked near the end of the interval, the accuracy 
deteriorates with increasing rriQ. The same, in principle, applies to the high exci- 
tations of light hadrons. A more appropriate strategy for heavy quarks is to start 
with a solution of the infinite-mass (static) equation. This has been done analytically 
1^, H, and full consistency^ with the OPE was demonstrated. 

However, our practical interest lies in the properties of heavy hadrons with rriQ 
lying in the intermediate domain, specifically for rriQ one order of magnitude larger 
than j3. The convergence of the l/mq expansion in this case is too difficult to quantify 
analytically. This is just the situation where the numerical computations are best 
employed. 

Therefore, an important element of the analysis is to check the accuracy of the 
numerical computations of both the heavy hadron masses and wavefunctions at dif- 
ferent values of mq. To this end, we compute the masses and certain moments for the 
ground and first excited states, and compare them to the analytic l/mq expansion. In 
general, the terms in the 1 / mq expansion depend on a number of expectation values 
in the static limit, like the kinetic one fil = (Q{iDyQ), etc. However, one can show 



Tl| that the parameters appearing here through high order in 1/mq can be expressed 



in terms of just the asymptotic value A = Mhq — mq and the corresponding decay 
constant. These quantities are the ones most accessible to numerical evaluation; in 
particular, A is expected to be the most accurately determined quantity. 

Our main computations refer to the case of m^^ = 0.56/3, as chosen in Ref. 0. 
It corresponds (see Sec. ^ to a mass of the strange quark in QCD. The choice of a 
noticeable light quark mass may be motivated by an attempt to mimic the effect of the 
transverse gluons absent in D = 2, which in a certain respect supply some effective 
mass to the light quark. Clearly, this can be only a rather crude approximation, 
since the bare quark mass breaks chiral invariance. One can suppose, nevertheless, 
that this side effect is not too important for our purposes. The chiral symmetry is 
spontaneously broken anyway, and the presence of a massless versus a massive pion 
does not seem to be essential for the range of problems we address here. On the 
other hand, the effect of the transverse degrees of freedom is known to soften the 



X ^ 1 singularity of the heavy quark distribution function |T^, |16|, 0, similar to the 
impact of the light quark mass in the 't Hooft model. The behavior of the distribution 
function affects the inclusive decays of the heavy quarks in an essential way. 

We also present some results for msp = 0.26/?, partly to explore light quark depen- 
dences of matrix elements and partly to investigate the beginnings of failure of the 
numerical solutions as mgp —>■ 0. The number of Multhopp modes used is 500; we 
considered smaller A^ as well to study this dependence, but since the behavior was 
found to be stable, we do not dwell on it further here. 



^In the case that the fermions / created by the weak current have m/ = 0, this agreement was 
shown up to and including ©(l/mg) terms in the weak decay width in B, while terms up to and 
including 0{rnj/mQ) were shown to coincide with those in the OPE in In the current work we 
take TO/ = 0. 
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The masses of heavy hadrons obey [18 
where 

III = (QitDfQ), pI = -^(q(D-E)Q 

pL = -\{tT{Q{zD)'Q{x),QizD)'Q{0)})^^^, (7) 

and these expectation values refer to the infinite mass hmit. In the p^^ expression, q is 
the momentum variable conjugate to x, and diagonal transitions within the correlator 
have been removed. 

Since A in QCD traditionally denotes the mass difference between a ground-state 
pseudoscalar meson and its corresponding heavy quark in the large tuq limit [as it is 
defined in Eq. (|)], and we need it for a number of the excited states Hq ^ as well, we 
assign the notation 

A^^^^e^'^), (8) 

and use e and A throughout the paper on equal footing. Equations ©-(ll) hold 
for each state Hq^ with n = 0, 1, . . ., so that an implicit superscript (n) is to be 
understood in these expressions. 

According to Ref. fill], the following relations hold in the 't Hooft model: 



Ptt — n ' Pd — 5 Pnn 



8A(A^-m^+/5') + 3/3^^2 



36 

(9) 

Here F is the scaled decay constant in the heavy quark limit, i.e., 

roo fTT-Q 1 / t \ 

= / dt^n{t) = lim / dt^=0„ 1 = lim Cn^/W^, 



(10) 

where the superscript is suppressed if there is no ambiguity, 

dxipnix) , (11) 

and the exact relation between c„ and the decay constant of the nth excitation is 
given in Eq. (|68D . In the heavy quark limit one has 

A = mQ(l-x), pl = ml{{x')-{xf) , (12) 

but there are ©(l/mg) corrections to these relations. For further applications to the 
decay widths we also consider the scalar expectation value P 

QQ) = ^C-). (13) 



2M„^ M„^ \x 



Then the following expansions hold: 



mn = A \ — — 5 V U 

^ ^ ' tuq 24ml 



/1\ _ ^ Pd-P^ 



TTTV 



Mhq \x/ 2ml 




0\^ \ . (14) 



We note that values of /x^, pf,, or p^^ determined from the expansions Eqs. (0) 
suffer degraded numerical accuracy compared to those taken directly from Eqs. (|^) 
since A and F are determined from more stable expansions (in particular, they do not 
depend upon close numerical cancellations). Therefore, we use Eqs. (|^) as primary 
information and relegate Eqs. (14) to numerical checks. Our method of determining 
A from the M^^ — mq expression, designed to minimize the influence of potentially 
large uncertainties at large mg, is described in Appendix |C[ 

Values of M^^ — mq and and the averages in Eqs. (|T^ as functions of mq from 
m = 0.56/3 to 50/3 are presented in Table |l| for both the ground and first excited 
states. Similar results for just the ground state with m = 0.26/? are presented in 
Table ^ Based upon the 10 data points presented in Table |1| for the ground state, 
one may fit to a polynomial in 1/mq, obtaining 

1 3 3'^ 

-iMnn-mn) = 1.317- 0.086— - 0.050-V + O 

/? ^ mq mf 

^ = 2.032 - 2.775— + O ( At 




^(1-x) = 1.316-2.491— + 3.. , 

P mq ml 





2 

m 



= 0.8074 





^C-) - 1 + 0.0994- 0.0444 + 0(4). (15) 
Mhq \x/ ml ml \ml) 

The corresponding expressions using the approach of Appendix y (neglecting the one 
for Mff^ — mq, which is used as input and hence is identical through 0{P'^/mq)) read 

= 2.035 - 2.816^ + (4 h 
mq \mlJ 




^(1-x) = 1.318 - 2.544^ + 4.5124 + of4 
P mq ml \ml 
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'^((x')-{xf) = 0.8078 - 3.996 A + o f ^ 
/32 V\ / ; rriQ \mi 



, , 1 + 0.096^- 0.066^ + 0^ . (16) 
Mhq \x/ m^Q m-^Q \m}^) 

This agreement between the two approaches is quite excellent and is exhibited in 
Figs. |l|-§ for Mhq — niQ and the quantities in Eqs. ([T^); in general, the exact results 



are presented as points on a solid hne, while each fit using Eqs. (0) is presented as 
a dashed line. In Fig. ^ the analogous expression Mh^ — tuq for the m = 0.56(3 first 
excited state is presented, while Fig. uses the same methods and values from Table ^ 
to present MnQ—mq for the m = 0.26/3 ground state. In Fig. |I] and especially in Fig. ^, 
the quality of numerical results is seen (as expected) to begin breaking down at large 
rriQ and small m, since A^ = 500 is fixed. We conclude that the numerical routine we 
rely upon is sufficiently accurate for = 500 up to niq ^ (25 ^ 30)/3. The critical 
value of TJiQ also depends, however, on the meson's light quark mass, decreasing for 
small m. This is expected since at small m the sharpness of the wavefunction as 
X — i> 1 becomes stronger, and more Multhopp functions are required to approximate 
it: Each Multhopp function vanishes like \/l — x. Likewise, the required N increases 
for the excited states. Still, one can check that it is possible to go as high as mg = 15(3 
even for m as small as 0.1(3. 

It turns out that a numerically significant cancellation occurs in the value of ii^—(3'^ 
in l/mg corrections and, in particular, at the l/m?Q level between p|, and p^^, for the 
ground state just around our primary value m = 0.56/3. Such a numerical suppression 
of the power corrections is accidental and does not occur for the excited states, nor 
for m = 0.26/3. 

Let us note that the expectation value of the light-quark scalar density in the 
heavy meson turns out very close to unity for the ground state, which may be seen 
by taking tuq ^ m and x ^ 1 — x in Eq. (|I^) and referring to Table 0; this is 
a characteristic feature of a nonrelativistic (with respect to the light quark) bound- 
state system. It implies an almost simple additive dependence of A on the light quark 
mass m, 

A ^ A|„=o + m , (17) 

and indeed one can verify this feature by comparing MnQ—mq values between Table |l| 
(m = 0.56/3) and Table | (m = 0.26/3). While this pattern is expected when the 
spectator quark is heavy, it a priori needs not hold when it is light. This supports 
the naive expectation that the chiral symmetry breaking may lead to a description 
in some aspects resembling the nonrelativistic constituent quark model. The above 
expectation value, however, decreases for the excited states, as expected from such a 
picture. 

Drawing semi-quantitative conclusions for QCD requires a translation rule be- 
tween the mass parameters in the two theories, that is, an estimate of the value of 
(3 in GeV. Different dimensionful quantities can be taken as the yardstick; since the 
theories are not identical, this translation rule must be introduced with some care. 
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As follows from the heavy quark sum rules, the physics of duality in the decay widths 
of heavy flavors crucially depends on the properties of the lowest excited heavy-quark 
states, in particular the P-wave excitations with opposite parity to the ground-state 
multiplet. It will become evident from the next section that they are of primary 
importance for the 1 / mq expansion of static properties as well. Therefore, we choose 
the mass difference between the lowest parity-even (P-wave) state and the parity-odd 
ground-state meson to gauge the translation between the mass scales. 

In the 't Hooft model the mass difference ei— eo for light spectators amounts to 
about 1.3(3. Real charm spectroscopy suggests that the first P-wave excitations are 
between 400 and 500 MeV above the ground state. Taking the larger value for sake 
of illustration, we arrive at the estimate 

p ^ 400 MeV, (18) 



which is adopted in our analysis. This falls rather close to the estimate of Ref. [p!9|] , 
which relied on a quite different type of effects in the light-quark systems. 

Assuming a value for the "bare" b quark mass in QCD (normalized at the appro- 
priate scale > nib) of about 4 GeV, we conclude that mesons with quarks of masses 

mb^lO/?, ^ (2.5^3. 5)/3 , (19) 

represent in the 't Hooft model the actual beauty and charm mesons. 

The value of A ~ 1.3/3 ~ 500 MeV seems to be in a reasonable correspondence 
with the size of this difference in QCD when it is normalized at a low hadronic scale, 
A(l GeV) ^ (600 ± 60) MeV 0. 



It should be noted, however, that the kinetic expectation value in the 't Hooft 
model turns out to be rather small, fx^ ~ 0.8/3^ ~ 0.12 GeV^. This is not surprising, 
since the chromomagnetic field is absent in two dimensions, while it was shown 



2T| , |22[] to be crucial in the real case. Indeed, the comparison is better justified for the 
difference /i^— /i^. in actual QCD versus the value of /i^ in the 't Hooft model. These 



questions were discussed in detail in Ref. [P3, and can be easily understood using 



the sum rule representation. Due to the absence of spin in two dimensions, there is 
no difference between the would-be spin-1/2 and spin-3/2 light degrees of freedom. 
In particular, labeling the "oscillator strengths" r defined in the next section by spin 
rather than excitation number, ti/2 = T3/2 and ei/2 = ^3/2 effectively hold. Then the 
sum 

f^l = 3E 4ki/2p + el\T,/,\' el\n/,\' , (20) 

n n n 

and the latter sum is just the general expression for fi'^—fiQ: 

fxl-l^l = 9j2'l\nM'- (21) 

n 

Accepting such an identification suggested in Ref. |^ and the estimate n'^ — jj^Q ~ 
(0.15 ± 0.1) GeV^, we again observe a reasonable agreement with the findings of the 
't Hooft model. 
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4 Duality in the SV Sum Rules 



A useful theoretical limit — the so-called small velocity (SV) regime — was suggested in 



the mid 80's |]2J] as a theoretical tool for studying semileptonic heavy quark decays. 
This refers to kinematics where both b and c quarks are heavy, but the energy release 
is limited, so that the velocity of the final charm hadron is small. At large energy 
release the OPE for the width must converge rapidly to the actual hadronic width. 
Still, at fixed energy release the deviations, although l/mg suppressed, are present 
regardless of the absolute values of masses. 

In the SV regime the semileptonic decays proceed either to the ground-state charm 
final state, D or D* (the semi-elastic transitions), or to excited "P-wave" states of 
the opposite parity. Other decays are suppressed by higher powers of velocity, or by 
heavy quark masses. 

The equality of the sum of partial decay widths and its OPE expansion is achieved 
through the sum rules that relate the sums of the P-wave transition probabilities, 
weighted with powers of the excitation energies, to the static characteristics of the 
decaying heavy hadron. The onset of convergence of the OPE expansion for the 
widths is then directly related to the pattern of saturation of the sum rules by the 
lowest excitations. If higher states contribute significantly, they delay the onset of 
duality, while their absence leads to a tight quark-hadron duality after the first P- 
wave channel is open. 

Knowledge of degree of saturation of the heavy quark sum rules is also important 
for another reason: It determines the hadronic scale above which one can apply the 
perturbative treatment to compute corrections or account for evolution of the effective 
operators. The lower this scale, the more predictive in turn is the treatment of the 
nonperturbative effects in the OPE. 

A recent review of the SV sum rules can be found in Ref. |^ (the perturbative 
aspects are discussed in more detail in Ref. [^). For most practical purposes ad- 
dressed here, one can consider the perturbative effects to be absent in the 't Hooft 
model. In particular, the heavy quark parameters do not depend perturbatively on 
the normalization point, and there is no need in the explicit ultraviolet cutoff to 
introduce a normalization point. The sum rules we address are 

pI-1 = E ^nk, (22) 
^Ak = E(en-efe)r„\, (23) 

{f^l), = E(^n-6,)Vi, (24) 

n 

{pI), = m^n-eufrl,. (25) 

n 

Here k and n denote excitation indices for the initial and final states, respectively (in 
practice only transitions from the ground state are interesting, so we limit ourselves 
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to A; = 0; in this case the index k is omitted). The so-called "oscillator strengths" 
T parameterize the transition amplitudes into the opposite-parity states in the SV 
limit, 



2mr 



(26) 



where v is the velocity of the final state hadron. In the diagonal transition pi is the 
slope of the Isgur-Wise (IW) function of state \k): 



1 



2mQ 



(27) 



The expressions for Tnk and pk in terms of the light-cone wavefunctions are 



oo 

dt^n{t)t-'^k{t) 

at 



/■I d 

lim / dxipn{x) {l-x)—ipk{x), 



Pk 



dt 



lim / dx 

rriQ^oo Jq 



dx 

'^~^^dx ~ 2 



{2i 



The finite-mg corrections to the x-integral forms turn out to be rather significant, 
leading to significant problems in precision numerical studies. To avoid this problem 
we use the analytic expression for the inelastic amplitudes obtained in Ref. [HI 



'Tnk 



n 



d 
'dt 



2{en-ek) 



p(n)p{k) 



■l-(-l) 



(29) 



where F^") are the asymptotic values of the decay constants c„ scaled up by the 
factor ^mg, as in Eq. (p!oD. The constants F*^"^ are computed as the values of 
Cn yJf^Q at mg = 15/? (scc Table 0) augmented by the 1/mg corrections detailed in 
the first of Eqs. (p!4D , while values of e„ are computed using the procedure described 
in Appendix |C[ 

The results of the computations for the case mg = 15/3, m = 0.56/3 are presented 
in Table ^. Our central result is a surprisingly good saturation of the sum rules: The 
first {n = l) excitation generates 99.4% of p^, 98.5% of A, 96% of p^, and even 91% 



of pj-). The rest is almost completely saturated by the second P-wave state 



n ■ 



3), 



where the cumulative values for the same quantities read 99.92%, 99.73%, 99.1%, and 
96.7%, respectively. 

In terms of absolute numbers, the sum rules Eqs. (|22|)-(p5D would give p^ — 1/4 = 
0.529, A = 1.278/?, pi = 0.782/3^, and p% = 0.99/33, ^j^g ^^^^ ^Yiich gives = 
1.99v^, in fine agreement with the values obtained from the values obtained in the 
previous section via the methods described in Appendix |C[ The few-percent discrep- 
ancy corresponds to the accuracy in determinations of squared decay constants. 

The level of saturation by the lowest open channels is extraordinary. The explicit 
reason for such a perfect saturation of the sum rules involving even rather high. 
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~ e'^ powers of the excitation energy can be read off Eq. (^) — r's are inversely 
proportional to the third power of the excitation energy. With the asymptotics e„ ~ 
~ n / , the first excitation energy ei — eo is notably smaller than the next 
one €3— eo including three energy gaps. The general peculiarity of the 't Hooft model 
leading to such a saturation is not understood completely. 

With this pattern of saturation of the SV sum rules for the ground-state meson, 
one expects an early onset of the accurate duality for the inclusive widths in the 6 — *• c 
transitions, only slightly above the threshold of the first excitation. Demonstrating 
this result through direct evaluation of the decay widths is one of the purposes of the 
next section. 



5 Local Duality in the Decay Widths 

The semileptonic widths described in this work were considered in detail in Ref. 
Here we recapitulate a few basic points. The weak decay Lagrangian is 

^weak = {^l^y) ■ (30) 

In terms of the previous notation, Q ^ b, q c (or, later in this section, u), and 
Hq — i> B. The key property of all D = 2 vectorlike currents is that for me = m,y = 0, 
the invariant mass of the lepton pair is always zero. For all computational purposes 
decays into this massless lepton pair are equivalent to decays into a single massless 
pseudoscalar particle weakly coupled to quarks according to 

£weak = ^c^^be^'d^cj). (31) 

V ZTT 

Several arguments favor our choice of a vectorlike weak decay interaction in the 
't Hooft model. One is of course the simplicity of Eq. (0). Another is that for 
q^ = some difficult problems of renormalization are absent, as we now discuss. The 
central problem in applying the OPE in practice is disentangling perturbative and 
nonperturbative effects. More precisely, this refers to the separation of short-distance 
effects attributed to the coefficient functions from long-distance effects residing in the 
matrix elements of the effective heavy-quark operators. 

The perturbative corrections, for example those that renormalize the weak quark 
current, are generally rather nontrivial, even in the 't Hooft model. However, accord- 
ing to the nonrenormalization theorem of Ref. such vertex corrections are absent 
from the decays with q^ = 0. This allows one to isolate the problem of renormalization 
of the underlying current from the question of interest in our study: possible devia- 
tions of the full decay widths due to the presence of thresholds in the production of 
the hadronic resonances. 

In reality, from the OPE viewpoint some short- distance corrections still remain 
even in this special kinematic region due to the high-momentum tails in the me- 
son wavef unctions. These tails come from the hard gluon exchanges between the 
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constituents. In principle, these "hard" components can also be separated from the 
"soft" bound-state dynamics explicitly. However, in practice this is not necessary: 
These effects are completely contained in the meson wavef unctions. 

Another advantage of vectorlike currents is apparent when one notes that the 
D = 2 three-body "semileptonic" phase space diverges logarithmically for massless 
leptons. Explicitly, for the decay M — »• m+mi+mi (equal lepton masses are assumed 
to render the expressions simpler), the three-body phase space turns out to be 



$3(M; m,mi,m(,) 



47r3(M-m)^(M+m)2-4mf 

(M+mf [(M-m)2-4m|] 
(M-m)2 [(M+m)2-4mf] 



xK 



(32) 



where K is the complete elliptic integral of the first kind. As mi —>■ m, one regains 
Eqs. (4.1)-(4.3) of 0], while as nii 0, the argument of the elliptic integral goes to 
unity, and K{l — e) ln(8/e)/2. This is a manifestation of the logarithmic infrared 
divergence of the massless scalar Green function at large distance in D = 2. A detailed 
calculation shows that the vector nature of the weak coupling regularizes the phase 
space integral, preventing the partonic rate from diverging in the limit of massless 
leptons. Furthermore, as discussed in the Introduction, this also removes the l/\p\ 
singularity in the threshold behavior for hadronic two-body decays. 

As a final advantage of vectorlike currents and the special kinematic point q^ = 0, 
note that at = the B — > D^''^^ transition amplitudes are directly expressed in 
terms of the overlap between the initial and the final wavefunctions: 



1 



5/^ 



2M 



{n\e 



B 



Jy\B) = -q^ ( <lxipn{x)ipB{,x) 
Jo 



(33) 



where Qz = —\p\ = 
D^^^ iu is given by 



-(M| 



■ M^)/2Mb, so that the partial decay width for B 



Ml- Ml 



Air 



B 



e{MB - M„) 



(34) 



The threshold suppression mentioned above is manifested in the explicit factor {Mi- 
Mi): The reciprocal of this factor in the phase space is removed by q1 from the 
matrix element. It is also possible to derive this result directly using the methods of 
Ref. 10; note, however, that these expressions are much simpler than those of Ref. ||^, 
because the vectorlike current with massless leptons restricts to 0. The sum of 
these widths over all open channels is to be compared to the OPE prediction. The 
remarkable speed of saturation in n, anticipated in the last section, is illustrated for 
one sample case in Table ^. 

Turning to the OPE, we mention one more problem associated with an accurate 
understanding of local duality violation. Apart from the purely theoretical aspect 
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that OPE power series are generally only asymptotic and, thus have a formally zero 
radius of convergence in l/mq, one normally has additional practical limitations. 
Only a limited number of the terms, as well as the associated expectation values, 
are usually known, which places additional theoretical uncertainties that dominate 
in practice at sufficiently large rriQ. 

This feature can be naturally incorporated in the analysis of our concrete model. 
We account completely only through terms that scale like l/rriQ, the highest or- 
der that emerges from the OPE free from the four-fermion operators 0. The rest, 
although calculable in principle term-by-term in the 't Hooft model, are taken to 
represent the OPE "tails" discarded by the unavoidable truncation. 

Using the sum rules of the 't Hooft model, Ref. P] established the following exact 
representation for the total decay width: 



rB = — - 



An 



2 2 



M, 



B 



— X 

X 



E 

M„>Mb 



(35) 



where r„ at M„ > Mb are understood as given by Eq. (|3^ ) without the explicit 9- 
function singling out the open channels; such are therefore all negative. On the 
other hand, the OPE yields the result 



nib 



nib 



dx 



Mb Jo x 




(36) 



with M generically denoting the OPE expansion parameter; we do not specify here 
if it is nib or mb — rric, or some other combination. It was shown in Ref. that the 
r„ term in Eq. ( |55D is dual to the order term in Eq. (P^D; however we do not use this 
here and rather treat the latter as an intrinsic uncertainty in the "practical" version 
of the OPE. 

Thus, our strategy is to compare the exact width 



B 



E 



Ml 



-Ml 



4^ Klt^Ms 



(lxipn{x)(fB{x) 



(37) 



to 



r 



OPE 



G^ 

An 



mt 



-m: 



nib 



nib 



1 dx 



X 



iflix) . 



(3^ 



The expectation value 



Mb Jo 

above can either be evaluated numerically, or in the 



spirit of the OPE, computed in the form of a l/m^ expansion, the last of Eqs. (0). It 
turns out that the expansion converges very rapidly to the exact result, so that this 
does not significantly affect the observed pattern of local duality at the quantitative 
level. The Born-term partonic rate is simply given by Pope with this expectation 
value set to unity. 



P 



mt 



-mz. 



nib 



(39) 



16 



The main practical interest of these calculations lies in the b ^ c width with 
its limited energy release E^. In general, can be small either if rrif, is not large 
enough, or even at large nib if E^ = rrih — rric (or mb — rric — \fq^ if is nonzero) 
is insufficient due to a significant c quark mass. The latter case falls into the SV 
category, and the violations of duality are suppressed here even at the maximal 
by heavy quark symmetry, as was pointed out in the mid-80's p4|. Therefore, one a 



priori expects a different pattern in the two cases. We try to separate the possible 
effects by considering different choices for and mc rather than by only taking them 
close to their realistic values. 

With these arguments in mind, one can expect to find significant effects of duality 
violation in the cases where l/m^ or 1/mfe effects are important. As suggested in 
Ref. [P3, in this case it is advantageous to fix close to its actual value, and vary 



mc from near m?, down to smaller values, changing in this way the energy release. At 
one end of the interval the local duality is supported by the heavy quark symmetry 
with large quark masses and SV kinematics, while at another end it rests on the large 
energy release. 

We start from the SV case when mt is fixed and large and vrtc is large as well, 
varying the energy release by increasing rUc towards mf,. Since the violation of local 
duality is expected to be suppressed for all values of mc, high numerical accuracy 
is vital. We fix m^ = 15/? (~ 6 GeV), and vary vtic from 5/? up to m\,. The results 
are given in Table ^ and Fig. ^. We note that the difference between the two widths 
is so small that one must plot ln(rB/ropE~ 1) rather than the widths themselves. 
This is expected since the SV sum rules are very well saturated, as detailed in the 
previous section — the higher thresholds are then strongly suppressed numerically at 
finite energy release. But for vtlc approaching m?,, where they could be noticeable, 
the heavy quark symmetry works efficiently since both quarks are very heavy. In 
fact, the only prominent features on the plot occur when thresholds to the first few 
D states of opposite parity to the ground-state B meson are crossed, for example 
between mc = 13.5 and 14/5. The deviation is extremely small also for smaller mc 
where the c quark velocity is rather large — yet there the energy release is significant, 
and a large number of excited states (up to 18 at mc = 5/? ~ 2 GeV) are produced. 
Table |^ and Fig. ^ show analogous results for mfe = 10/3, m = 0.56/3. 

To render the duality violation more apparent, we consider (Table |]and Fig. 0) 
the same decay widths for a h quark with half the mass, m^ = 5/3 ~ 2 GeV. Even 
here the deviation is below per mill as soon as the first excitation can appear with 
sufficient phase space. The duality-violating component at last exhibits the proper 
oscillating behavior (note the decrease between mc = 3 and 3.5/? or 1 and 1.5/3), but 
this effect is too small to be extracted reliably at larger energy release where this 
property becomes an asymptotic rule. 

As follows from our computations, local duality is violated at a tiny level in the 
6 — s> c decays in the 't Hooft model whenever it is a priori meaningful to apply OPE. 
A possible reason behind this might be that for unidentified reasons the heavy quark 
symmetry works for the inclusive widths too effectively, down to relatively low masses 
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and velocities of order 1. This was conjectured in the early papers on the subject 



24|. 

Therefore, our final attempt in the quest for the sizeable duality violation in 
beauty is considering the (6^u)-type transitions, where the heavy quark symmetry 
per se does not constrain the individual transition form factors. We fix in our ex- 
pressions mc = m = 0.56/3 or 0.26/3 (but still keep the two quarks flavor-distinguished) 
and vary rrib from l/?^0.4GeV to 12/3^^4.5 GeV. The results are shown in Table ^ 
and Fig. |Tl], and Table |10| and Fig. respectively. Although the difference between 
the actual width and its OPE approximation is larger, it still is very small and ap- 
proaches a percent level for as low as 2(3^0.8 GeV. The total decay width is no 
longer saturated to such a high degree by transitions to the ground state, especially 
for larger rrih. Nevertheless, the duality is amazingly well satisfied when just the 
first few open channels are summed. Again, the only prominent features in the plots 
appear when crossing kinematic thresholds due to the lightest D mesons of opposite 
parity to the ground-state B. 

The extraordinary agreement between and Fqpe may be underscored by in- 
stead plotting (Fig. |13|, final column of Table the difference between F^ and the 
Born-term partonic rate F^ given in Eq. (|39|) . From an algebraic point of view, F^ 
and Fqpe differ generically at 0(1/M^), while F^ and F;, begin to differ already at 
0(1/M2). 

Thus, we find local duality between the actual semileptonic decay width and its 
OPE expansion to be very well satisfied in all cases. 

Before concluding this section, let us briefly address duality in the differential 
distribution F~^ dT/dE. In the heavy quark limit the shape of the final-state hadronic 
mass distribution follows the heavy quark distribution function in the decaying meson; 
for the 6— decays under consideration, this is the light-cone distribution function 
F{x). In decays with g^ = the recoil energy of the lepton pair E is directly related 
to the final state mass M^. 

^ Ml- Ml , ^ 

Since = 0, these decays are analogous to 6 ^ 57 in the Standard Model. In the 
large-mf, limit one has 

IdF ^ 2 /2E-^X 

FdE A V A y ^ ^ 

At finite mf, in a theory with narrow resonances the actual distribution is given by 
the comb of 5-functions with spacing in the argument of Eq. (^) of order A/m^. 
In order to obtain a continuous result, we adopt the simple ansatz of averaging 
over the peaks. Using Eq. (|40| ) to define the energy En of the nth state Mh, we 
integrate the (5-function for the nth state evenly over the energy range {En+En+i)/2 
to {En + En-i)/2, i.e., the midpoints between energy eigenvalues. Letting be the 
maximum number of kinematically allowed Mh values, we establish the endpoint bins 
by defining E_i = E^ax = Mb/ 2 and E^+i = ^min = 0. 
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We find that our numerical computations yield a distribution resembling the light- 
cone distribution function specifically, 

= A*'((l-y)A) = ^Ijm^A^^fl- 

Recalling that M5 = m;, + A+0(l/mfe) and combining Eqs. 

1 dr 20/. Mb-2E 




The two sides of this expression are plotted in Fig. using rriQ = 25/3 to represent 
the limit mg — *■ 00, while the actual distribution is considered at rrib = 10/3. The 
agreement is quite remarkable. The continuous distribution appears to pass approxi- 
mately through the midpoint of each bin; owing to the near-equal spacing of 't Hooft 
model eigenvalues in M^, Eq. ( ^O]) shows that these bin midpoints are very close to 
the values En themselves. 

It is also interesting to consider integration over a range of E. In particular, define 
^{1~2E/Mb) as the cumulative fractional width from maximum energy Mb/2 down 
to the given E; then $(0) = and $(1) = 1. While the exact result for ^r/r 
amounts to an integration of the 5-function differential widths renormalized so that 
the cumulative result approaches unity, the integral of the continuous distribution 
gives 

2 fv . , / Mh^ ' 



$(?/)= lim / d;Z¥?Ml ^z\. (44) 

mg-^oo rnq Jo \ J 

These two curves are presented in Fig. |l^. Two features particularly stand out in 
this plot. First, even for rrih as large as 10/3~4GeV, the overwhelming part of the 
decay probability falls into the transitions to at most four lowest states. Second, 
the continuous curve seems to provide a nearly optimal description possible for the 
step-like exact distribution. The point-to-point deviation for all plotted values with 
1 — 2E/Mb>0.04 does not exceed half of the contribution of the nearest threshold. 



6 Duality in the Vacuum Current Correlator 

In this section we briefly illustrate the onset of duality for the absorptive part of the 
vector current correlator with light quarks, of the type that determines the normalized 
cross section R{e^e~ — > hadrons) as a function of energy. In the context of the heavy 
quark decays this is relevant in nonleptonic decay widths in two kinds of processes: 
in spectator-independent decays, where R{q^) determines the weight with which the 
semileptonic width at given must be integrated over (see Ref. 0), and in the 
effects of WA decays. 

In either case, at 00 the cross section appears as a comb-like collection of 

5-functions: 

R{q') = ll^lK<f-Ml) ; c„ = fdxipn{x) . (45) 
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The above expression for the residues refers to the case where a vector current is 



considered. We suppress here the factor of yNc/n relating c„ to /„ [Eq. (p8|)]. We 
also assume in what follows that the light quark masses are equal, mu = md, and are 
0{f3) or less, in order to reach asymptotic more quickly. 

In the extreme situation of infinitely narrow resonances one cannot, of course, 
discuss a point-to-point equality of the cross section R{q^) with its OPE in the form 
of expansion. A meaningful comparison is possible if each resonant peak is 
somehow averaged over an interval no smaller than the distance between adjacent 
peaks, the latter being approximately given by Ag^ ^ 7r^/5^ [|l|. It is worth recalling 
that R{q^) is proportional to rn? jq^^ so one must consider nonvanishing masses for 
the vector current, and address the OPE terms formally suppressed by m? /q^. 

This question was first addressed in the context of nonleptonic decays in Ref. 
using the numerical approach. Duality for the average cross section in the same man- 
ner as above, i.e., using sum rules derived from the 't Hooft equation and analytically 
matching terms in the l/mq expansion, was obtained in Ref. [^. Yet establishing 
the asymptotics per se cannot tell us beforehand how early one can expect the onset 
of duality. Here we study this question numerically, in the domain of intermediate 

The concrete amount of the deviation between R{q^) and i?°^^(g^) in the case 
of direct resonances may depend in an essential way on the chosen smearing pro- 
cedure. Interested in the qualitative features only, we choose a rather simplified, 
crude method: We spread the integral of Riq"^) evenly over the interval between the 
successive resonances. More precisely, we put 



Ri^') = 772 ^77^ Z"""^' dq'R{q') = „ . "^'^ . , (46) 

for M2n-i < q^ < M2n+i, with M^^ = 4m^, the partonic pair production threshold. 
Here we use the fact that c„ vanish for odd n when = rrid = m. This smearing 
is very similar to that described for the differential width in the last section, except 
that averaging is performed in rather than E. The free quark loop R{q^)i which 
is of course the leading term of the OPE, is given by 

„ , o, 2m? 1 

q^ y/l-4mVg2 

Table [TT] and Fig. [T^ show the results for our reference case m = 0.56/3. The agreement 
of the average hadronic cross section with the parton-computed probability again 
turns out to be very good. Apparently, this can be related to two facts: the heavy 
suppression of power corrections to -R(g^) in the OPE (see Eqs. (34)-(35) in |^), and 
an early onset of the asymptotics in the spectrum, 

M„V-M„2 ^ (48) 

which even at n = is satisfied to about 15%. 
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7 Discussion and Summary 



The main motivation behind the present study has been to assess the magnitude of 
local duality violations in the inclusive semileptonic decays of beauty particles. We 
considered this question using the 't Hooft model as a toy theory in which all relevant 
decay amplitudes can be evaluated numerically. The 't Hooft model, while retaining 
certain key features of full D = 4 QCD that shape the spectrum of hadrons (quark 
confinement, chiral symmetry breaking), still differs from D = 4 in many respects. 
Yet using it as a lab for exploration carries an important advantage — it allows no 
"wiggle room" for interpretation of the results. There are no ad hoc parameters to 
choose or adjust, and as soon as the underlying weak decay Lagrangian is fixed, the 
numerical results are unambiguous and must be accepted at face value. This posi- 
tively distinguishes this approach from various models where often the conclusions, 
even qualitatively, depend on the arbitrary choice of parameters according to one's 
preferences. The question of a particular model being compatible with the general 
dynamical properties of QCD underlying the OPE approach, often quite problematic 
in simplified quark models, does not arise for the 't Hooft model. 

Although the simplest illustration of the asymptotic nature of the decay width 
l/mg expansion and related violations of local duality ||^ follows just from the ex- 
istence of hadronic thresholds (see, e.g., violation of local duality is a more 
universal phenomenon that is not directly related to existence of hadronic resonances 
nor even confinement itself. This has been illustrated in Ref. [Q by the example of soft 
instanton effects that do not lead, at least at small density, to quark confinement — 
but do indeed generate computable oscillating duality- violating contributions to the 
total decay rates. 

Nevertheless, there is a widespread opinion that decays with manifest resonance 
structure in the final state are most difficult for — if compatible at all with — the stan- 
dard OPE. Even the possibility that the OPE does not fully apply in the case of 
"hard" confinement has been occasionally voiced in the literature. The analytic 
studies performed in Refs. 0, which explicitly demonstrate in the 't Hooft 

model the applicability of the OPE to the total widths, should help to allay such 
conceptual concerns. Nevertheless, the intuition remains that resonance dominance 
is not "favorable" for the OPE, and problems might show up, for instance, through 
a delayed numerical onset of duality, in that the approximate equality of the OPE 
predictions and the actual decay widths may set in only after a significant number 
of thresholds has been passed. To address such issues, the 't-Hooft model seems to 
represent the most certain testing ground for local duality in the domain of decays 
of moderately heavy quarks. 

Contrary to naive expectations, we found surprisingly accurate duality between 
the (truncated) OPE series for Fgi and the actual decay widths. The deviations are 
suppressed to a very high degree almost immediately after the threshold for the first 
excited final state hadron is passed. No suspected delay in the onset of duahty was 
found. 
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The key property that governs the onset of the l/mg expansion for the semilep- 
tonic widths is the pattern of saturation of the heavy quark sum rules. We examined 
a particular class, the SV sum rules in the heavy quark limit, that has the most 
transparent quantum mechanical meaning. We found them saturated to an amazing 
degree by the very first excitation. The contribution of the remaining, higher states 
to the slope of the IW function. A, and /i^ does not exceed a few percent. Even in the 
Darwin operator sum rule, the first excitation accounts for 90% of the whole expec- 
tation value, despite the fast-growing weight, (e^— eo)^ of higher-order contributions. 
This peculiarity underlies the early onset of duality for the case when initial- and 
final-state quarks are both heavy. 

Some of the duality-violating features observed in these studies have natural ex- 
planations. At fixed energy release rriQ—mq the magnitude of the deviations is smaller 
if TTiQ, rriq are both large (as in b—>-c) than if they are both small. This is expected, 
since in the former case the heavy quark symmetry for the elastic amplitude addi- 
tionally enforces approximate duality even when no expansion in large energy release 
can be applied. 

It is interesting, however, that at fixed rrib the duality violation decreases rapidly 
as rric decreases, in full accord with the OPE where the higher order terms are gen- 
erally suppressed by powers of l/(mb— mc). This is clearly a dynamical feature that 
goes beyond heavy quark symmetry per se, the quality of which deteriorates as rric 
decreases. 

It is also instructive to note that including the calculated power-suppressed OPE 
terms significantly reduces the difference between the actual decay width and its 
purely partonic evaluation. Moreover, the seeds of oscillations inherent to duality 
violation (as functions of quark masses), can be seen. Since we adopted the trun- 
cated OPE expansion to mirror the existing implementation of the OPE in QCD, the 
deviations do not average to zero but rather oscillate around the (rapidly dissipating) 
contributions attributed to discarded higher-order terms. 

The numerical effects of duality violation we study turn out to be typically quite 
small. Partially this can be attributed to moderate size of the corresponding expec- 
tation values multiplying l/mq corrections in the OPE. Yet certainly not all power 
corrections in heavy quarks are suppressed in the model. It is well known from ordi- 
nary quantum mechanics that masses (eigenvalues) typically are much more robust 
against perturbations than wavefunctions themselves (or transition amplitudes). We 
observe a similar pattern in the 't Hooft model. For example, l/mq corrections to 
the meson decay constants turn out very significant even at the scale of the b quark 
mass. Apparently, the inclusive decay rates fall into the class of "robust" observables, 
although, as explained above, this was difficult to anticipate beforehand. 

We note here another "fragile" observable, the light-cone heavy quark distribution 
function, which can be measured in decays of the type b s'~f. In D = 2 the scaled 
spread ttiq ((x^) — (x)^^ of the x distribution approaches /i^ at large ttiq. Yet, as seen 
in Fig. ^, even at the b quark mass one would obtain from this distribution only about 
60% of the actual value of /x^, due to significant l/mq corrections. This caveat may 
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be important for existing analyses of the decay distributions in B decays, where such 
effects routinely are not included. 

We also briefly addressed the inclusive differential decay distributions in the ana- 
logues of b —>■ u iu or b —>■ sj decays. Generally, we find good agreement (at the 
scale corresponding to the physical b mass) with the parton-based prediction incor- 
porating effects of the "Fermi motion," and in particular for the partially integrated 
probability 

1 r^^^l dr 1 

= ^ / dM,2 1^ . (49) 
Tsi Jo dMfi 



This distribution, following Refs. |25, Eq], is examined in real B decays in the quest 



for \Vub\ 1231 • However, the point-to-point deviations are clearly still significant, for 
the decays to only the 4 or 5 lowest final states saturate the overwhelming fraction 
of the total decay probability. It is quite conceivable, though, that such deviations 
are less pronounced in actual QCD owing to the significant resonance widths and to 
a richer resonance structure. 

The vacuum current correlator also turns out to be especially robust; even neglect- 
ing all OPE corrections except the leading partonic contribution leads to excellent 
agreement with the hadronic result. 

Turning to the direct phenomenological conclusions that can be inferred from 
our studies, we see that, to the extent our findings can be transferred to real QCD, 
violation of local duality in the total semileptonic widths of B mesons is not an issue. 
The scale of duality violation lies far below the phenomenologically accessible limits, 
and cannot affect the credibility of \Vcb\ or \Vub\ extractions. 

In reality there are, of course, essential conceptual differences between the two 
theories, including those aspects that are expected to be essential for local duality 
(for a discussion, see Ref. |^). Although many seem to optimistically suggest that 
duality violation is more pronounced in the 't Hooft model than for actual heavy 
flavor hadrons, some differences may still work in the opposite direction. In D = 2 
there are no dynamical gluons, nor a chromomagnetic field that in D = 4 provides a 
significant scale of nonperturbative effects in heavy flavor hadrons. Likewise, there 
is no spin in D = 2, and no corresponding P-wave excitations of the light degrees of 
freedom (the so-called j = 3/2 states), which seem to play an important role in D = 4. 

Two-dimensional QCD neither has long perturbative "tails" of actual strong in- 
teractions suppressed weakly (by only powers of logs of the energy scale). In D = 2 
the perturbative corrections are generally power-suppressed, as follows from the di- 
mension of the gauge coupling. As discussed in Ref. it is conceivable that the 
characteristic mass scale for freezing out the transverse gluonic degrees of freedom 
is higher than in the "valence" quark channels. This would imply a possibly higher 
scale for onset of duality in a<j/vr corrections to various observables. 

Regardless of these differences, we conclude that presence of resonance structure 
per se is not an obstacle for fine local quark-hadron duality tested in the context of 
the OPE. As we see in the 't Hooft model, resonances themselves do not seem to 
demand a larger duality interval. As soon as the mass scale of the states saturating 
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the sum rules in a particular channel (quark or hybrid) has been passed, the decay 
width can be well approximated numerically by the expansion stemming from the 
OPE. 

The ground states of heavy mesons in the 't Hooft model exhibit relatively small 
expectation values of nonperturbative operators (/x^, p%, but not A) compared to 
real QCD, if our identification f3 ~ 400 MeV is adopted. This may be regarded as 
a reason for small duality violation for Fgi in the model. However, even if we scale 
/3 up to 700-800 MeV to make up for smallness of the nonperturbative OPE effects, 
the duality violation is still very small, and superficially rather insignificant even in 
charm. 

We note, however, that the specific choice Eq. (|30|) of the weak interaction effec- 
tively requires decays to occur only at q'^ = 0, and therefore the effects of four-fermion 
operators of the type {QTq){qTQ) are totally absent, at least in the lowest orders of 
perturbation theory (cf. Ref. 0, Sec. III.B.3). As was suggested in Ref. it is 
conceivable that the apparent excess in rsi(-D) is simply related to a noticeable mag- 
nitude of the non- valence (nonfactorizable) expectation values {D {{cTs) (sFc)! D). If 
this conjecture is true, similar effects in Tsi{B — > X„£z/) are still suppressed but possi- 
bly detectable in future precision experiments. In the context of the present study, it 
suffices to say that this would be a legitimate OPE effect rather than a manifestation 
of a significant local duality violation in the strict sense. 
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A The BSW Improvement of the Multhopp Tech- 
nique 

The Brower-Spence-Weis [|I^ (BSW) improvement of the Multhopp technique avoids 
the need for evaluating the wavefunction at a discrete set of points called "Multhopp 
angles," thus improving the behavior of the solutions in the endpoint regions, as 
described in Sec. ^. Here we exhibit the expressions used by BSW, correcting along 
the way some minor typographical errors in their work. 

Starting with the 't Hooft equation (0) with bare quark masses mi and m2, one 
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converts the kinematic variables x, y to angular variables: 



X 



1 + COS6' 



1 + cos^' 



(50) 



in terms of which the 't Hooft equation reads 



M2 



+ 



l+cos6' 1 — cos6' 
JO (cos fc* — COS 6'')"' 



Expanding 



Vp{0) = E 4"^sinn^, 

n=l 

and using the continuous inversion identity (contrast with Eq. (A7) of Ref. 0) 



d^ sinm^ sinn^ 



i,^mn ^m,,—n) 



one obtains the infinite-dimensional eigenvector system 



nm m i 



where 



4 r 

+- / d^ 

TT Jo 

-'-f r 

IT Jo 



+ 



rrin 



sinnO sin m6', 



l + cos6' 1— cos6' 

/•TT 

dO sin n9 / d^' sin 9' sin m^' P 



l^cos 6' — cos 6*')^ 



Both of these integrals can be evaluated, with the result 
{Ho)nm = 4min(n,m)[(-ir+"m? + "^2" 



m 



m—lj n+m — 1 



1 + (-1)"^ 



(51) 
(52) 

(53) 
(54) 

(55) 
(56) 

(57) 
(58) 



where Vno = Vom = for m,n>0. This recursive form for Vnm is most convenient for 
numerical calculations; however, one may also write the closed-form solution. 



4m 



1 + 



\n+m 



1 — n—m 



1 — In— ml 



(59) 



Note that the "potential" V in Eq. ( |55D is real but not symmetric, owing to the extra 
sin^' in Eq. (|56|); therefore, the "Hamiltonian" Hq + V is not Hermitian, and the 
eigenvectors a^^^ are not orthogonal. This is a direct result of converting the exact 
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wavefunctions, which are eigenfunctions of a Hermitian Hamiltonian when written in 
terms of the variable x (and therefore orthogonal in x), into orthogonal functions of 
the variable 6. This transformation is nonunitary because the number of modes used 
is not infinite; therefore, the overlap of different eigenvector solutions should be small 
when a large number of modes are used. Indeed, this turns out to be empirically true; 
nevertheless, we take the further step of orthogonalizing the numerical eigenvector 
solutions recursively by means of the standard Gram-Schmidt procedure, i.e.. 



A) 



ip)\ _ 



V'orth 



(P) 



(P) 



V'orth 



(60) 



For = 500 modes, this typically changes expectation values by one part in 10^. 

The expressions for these overlaps and other matrix elements in terms of the mode 
coefficients a„ are presented in Appendix p. 



B Matrix Elements 



A number of useful overlaps and other integrals are straightforward to evaluate in 
terms of the mode coefficients, using the expressions (^^. Solving them amounts 
to evaluating a number of trigonometric integrals. Such expressions are especially 
convenient since they permit a number of integrations that introduce no numerical 
uncertainties (except due to machine precision) beyond those of solving the original 
Multhopp-BSW eigenvector equation (|5^) . 

In particular, denote the pth eigenstate wavefunction presented in Eq. (|^ ) by 
if^^ and that for some other set of masses in the gth eigenstate by ip^^; the latter 
wavefunction then has an expansion like (|5^) with mode coefficients bl^\ Truncating 
after N modes, one then finds 



dxy^'^\x)y^^:Kx) 







N 

—2 ^ ma, 

m=l 



q 

N 
n=l 



l + i-l) 



[1 — (m — n)2] [1 — (m- 



-n] 



(61) 



Indeed, the normalization integral Jq dx ip{x)'^ = 1 is just the case a = b and p = q, in 
agreement with Eq. (A9) of Ref. ^j. 

Other useful expectation values include 

dx (x — 



N 



N 



na. 



m=l 



n=l 



l-(-l) 



m+n 



1 



[A-{m-nY] [4-(m+n)2]' 



(62) 
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X — - 



dx (x — 
V 2 



N 



N 



na 



ip) 



m=l 



n=l 



Vi^\x) 

l + f_l)m+n 



X [(l-(m-r2)^) (l-(m+n)^) (9-(m-n)^) {9-{m+n)'- 



-1 



(63) 



Note that the spread of the wavefunction may be computed about any convenient 
point in x, viz., 



{{ax+bf) - {{ax+b)y = a" 



so that the additive constants of —1/2 above are irrelevant. Also, 



Xi p 



N 



N 



dx — 

X 



(64) 



(65) 



■m=l n=l 



where 



(m+n)/2-l 



2 E 



j=\m-n\/2 



2j + l 



\m — n\ + l 



{m+n-l)/2-l 



m — n even; 



E 



\m — n\ m+n . ,,,„27 + l 



1 



One also finds 



1 — x 



\m—n\ m+n 



dx 



1 (\m — n\ 



( m+n \ 



m—n odd. (66) 



1—x 



TV 



= E E^lf^-^" 



(67) 



m=l n=l 

one finds J^n = +Imn for m — n even, and 



where, using the notation of Eq. (^ 

Jmn = -Imn for m-Tl odd. 

Finally, the decay constant of the pth excitation [cf. Eqs. (p!0D-(p])] is given by 



C Additional Relations Used in the Analysis 



(68) 



The numerical calculation of large-mg matrix elements with acceptable accuracy 
relies on achieving a balance between competing effects. 
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On one hand, Multhopp solutions to the 't Hooft equation with tuq ^ (3 tend 
to suffer degraded numerical accuracy since they are highly concentrated into the 
small kinematic region 1—x <^ 1. As discussed in Sec. 0, the endpoint regions x~0 
and 1 are where the Multhopp solutions — or more precisely, their derivatives — tend 
to break down. This effect is compounded when m <^ (3, since lighter quark masses 
force sharper endpoint behavior in the wavefunction. Although the BSW solution 
ameliorates this behavior, as rriQ is increased one eventually faces the problem of 
attempting to represent a function with only a very small region of support in x by 
a finite number of modes with support over the full range x G [0, 1]. In practice, we 
gauge the errors committed through such "lattice spacing" effects by computing a 
given quantity with = 500 and noting the amount by which its value shifts if one 
uses instead A^ = 100, and as expected, such errors become substantial (as much as a 
few percent) by the time one reaches mg >25P or m< 0.4/3. 

On the other hand, although numerical solutions with mq, m ~ 0{l3) have the 
highest numerical accuracy, they also have substantial O^I/tjiq), ©(l/mg), etc. cor- 
rections that are difficult to disentangle. 

We adopt an intermediate strategy of employing certain exact relations that hold 
for the 't Hooft solutions. To determine the relevant static expectation values, we 
solve the finite-mg heavy hadron mass expansion for A [Eq. (|^)]: 

Neglecting the order term and using the relations [Eq. (|^)] 



-r2 



P'lT ~ 2 ' Pd ~ ^ ' Ptttt 2g 



A{A -m'+P') + 3P'F' 



(70) 

we thus arrive at an equation cubic in A that depends on F^. We solve it at mq = 15/3. 

The asymptotic value of the scaled decay constant F^"'^ = ^Jmq Cn must also be 
evaluated at a finite value of ttiq, thus including l/mg-suppressed pieces. We account 
for them explicitly using the expansion |]ll| [the first of Eqs. (|l^] 

We likewise solve this equation for F^"^ at ttiq = 15/3. 

Turning to the analysis of the SV sum rules Eqs. (p^)-(|25|) in Sec. 4, we note that 
their rapid saturation demands an exceptionally high precision in evaluating both the 
oscillation strengths r in the r.h.s. and the expectation values in the l.h.s. Reaching 
such an accuracy through direct computation seems impossible. Therefore, we use a 
number of identities to get meaningful results. First, we employ the expression for 
Tnk in terms of e„, e^, and the corresponding decay constants: 
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Then we make use of the fact that the discussed sum rules, being completeness sums, 
are exact when summation includes all excitations (see Ref. UTTII). Therefore, one has 



Pk 4 



1 " 

2 _ _ \ ^ 2 

Pfc T ^ ''"^fe 
^ 1=1 



T E 



2_1 ^ ^^fc ) 
Pk 4 £=n+l 



2 

Afc 



-Afc -^(e^-efe) 



1 



ik 



2 _2 



D)k L 



oo 



oo 



{P' 



D)k £=n+l 



(73) 
(74) 
(75) 
(76) 



The sums on the r.h.s. can be accurately evaluated since the higher contributions fall 
off in magnitude very quickly. In practice, we truncate the sum at £ = 20. 

A similar approach was used to evaluate the duality- violating difference F^ — Tope 
as a function of mq. We use the exact relation p [Eqs 



), dH), (H)] 



^ G"^ ml — ml /"i dx , , , 

^B = -. \^ 

An Mb Jo X 

and therefore, 

rs-TopE / /"^ dx 



E 



Ml-Ml rl 



4^ M„>Ms Mb 



dxipn{x)(pB{x) 



, (77) 



OPE 



X 



E 



ml — m^ 



(lxipn{x)ipB{x] 



9{Mn-MB) 



(78) 

The summation runs over all final excited states kinematically forbidden in the decay. 
Once again, the sum converges rapidly and is dominated by the lowest couple of states. 
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Mng-rriQ 




mQ(l-x) 




"iQ /1\ 1 


Ground state 


n = 0) 


0.56 


1.21918 


0.7300 


0.280 


0.017 


8.34 X 10-2 


1.0 


1.24633 


0.9534 


0.432 


0.048 


4.59 X 10-2 


3.0 


1.28764 


1.4210 


0.791 


0.211 


9.04 X 10-3 


5.0 


1.29904 


1.6061 


0.944 


0.333 


3.58 X 10-=^ 


7.0 


1.30423 


1.7503 


1.029 


0.417 


1.88 X 10-3 


10.0 


1.30820 


1.7901 


1.102 
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9.33 X 10-^ 


15.0 


1.31131 


1.8633 


1.166 


0.582 


4.14 X 10-^ 


25.0 


1.31375 


1.9271 


1.222 


0.661 


1.46 X 10"^ 


35.0 


1.31475 


1.9560 


1.248 


0.700 


7.31 X 10"^ 


50.0 


1.31545 


1.9783 


1.268 


0.732 


3.53 X 10-5 
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2.82831 


0.0000 


0.280 


0.032 


-1.30 X 10-1 


1.0 


2.77888 


0.0922 


0.476 
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-9.62 X 10-2 


3.0 


2.66569 


0.4429 


1.094 


0.457 


-3.29 X 10-2 


5.0 


2.61977 


0.6427 


1.437 
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-1.60 X 10-2 


7.0 


2.59522 


0.7648 


1.649 


1.014 


-9.40 X 10-3 


10.0 


2.57436 
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1.848 


1.267 


-5.14 X 10-3 


15.0 
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0.9812 


2.033 


1.529 


-2.50 X 10-3 


25.0 


2.54088 


1.0765 


2.205 


1.797 


-9.68 X 10"^ 


35.0 


2.53382 


1.1213 


2.287 


1.933 


-5.10 X 10"^ 


50.0 


2.52833 


1.1566 


2.351 


2.046 


-2.56 X 10"^ 



Table 1: Matrix elements as functions of heavy quark mass mq and light quark mass 
m = 0.56/3, computed numerically via the BSW-improved Multhopp technique. All 
masses are in units of jS. 
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0.5067 
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3.61 


X 10^ 


-1 
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0.658 
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1.40 


X 10- 


-2 
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X 10- 


-3 
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1.6277 
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2.80 


X 10- 


-3 
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1.6999 
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1.38 


X 10- 


-3 


15.0 


1.02692 


1.7612 
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6.51 


X 10- 


-4 


25.0 


1.02946 


1.8137 


0.969 


0.591 


3.42 


X 10- 


-4 


35.0 


1.03003 


1.8369 


0.987 


0.665 


2.96 


X 10" 


-4 


50.0 


1.02992 


1.8543 


1.002 


0.811 


2.97 


X 10" 


-4 



Table 2: The same matrix elements as in Table |I]for the ground state and m = 0.26/9. 
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-8.3 X 10- 


-2 


1.00 


21.71902 






20.2 


-7.0 X 10- 


-2 


1.50 


22.18785 






14.0 


-5.6 X 10- 


-2 



Table 3: Ground-state matrix elements as functions of m, used to probe the expec- 
tation value of the light-quark scalar density {Hq \ qq \ Hq) /2Mhq ■ 



k 








1/4 




A 


/i 


2 

IT 




Pi 





1.318 






























1 


2.516 


7.25 X 


10- 


-1 


6.2 X 


10- 


3 


1.5 X 


10- 


-2 


3.6 X 


10- 


-2 


8.8 X 10- 


-2 


3 


3.989 


5.36 X 


10- 


-2 


7.8 X 


10- 


4 


2.7 X 


10- 


-3 


9.4 X 


10- 


-3 


3.3 X 10- 


-2 


5 


5.060 


1.73 X 


10- 


^2 


2.2 X 


10- 


4 


9.3 X 


10- 


-4 


4.1 X 


10- 


-3 


1.8 X 10- 


-2 


7 


5.949 


8.37 X 


10- 


-3 


8.7 X 


10- 


5 


4.3 X 


10- 


-4 


2.2 X 


10- 


-3 


1.1 X 10- 


-2 


9 


6.724 


4.92 X 


10- 


-3 


4.1 X 


10- 


5 


2.3 X 


10- 


-4 


1.2 X 


10- 


-3 


6.7 X 10- 


-3 


11 


7.421 


3.24 X 


10- 


-3 


2.1 X 


10- 


5 


1.3 X 


10- 


-4 


7.4 X 


10- 


-4 


4.3 X 10- 


-3 


13 


8.060 


2.29 X 


10- 


-3 


1.1 X 


10- 


5 


7.0 X 


10- 


-5 


4.4 X 


10" 


-4 


2.7 X 10- 


-3 


15 


8.653 


1.71 X 


10- 


-3 


5.4 X 


10- 


6 


3.6 X 


10- 


-5 


2.4 X 


10- 


-4 


1.5 X 10- 


-3 


17 


9.209 


1.32 X 


10- 


-3 


2.1 X 


10- 


6 


1.5 X 


10- 


-5 


1.0 X 


10- 


-4 


6.6 X 10" 


-4 


19 


9.735 


1.05 X 


10- 


-3 

























Table 4: Meson mass eigenvalues = Mj^^—niq and oscillator strengths r as func- 
tions of excitation number k, for m = 0.56/5. The value for a given nonperturbative 
matrix element for each k indicates the fractional amount remaining after saturating 
the corresponding sum rules Eqs. (p^)-(]25|) by states n with n < k. 
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n 













11.3082 


0.96426 


9.45 X 10- 


-1 


0.9448313 


1 


12.5744 


0.25996 


5.36 X 10- 


-2 


0.9984744 


2 


13.4495 


0.044347 


1.23 X 10- 


-3 


0.9997062 


3 


14.1780 


0.023927 


2.74 X 10" 


'4 


0.9999801 


4 


14.8095 


0.005462 


1.03 X 10" 


'5 


0.9999903 


5 


15.3810 


0.006617 


9.49 X 10" 


'6 


0.9999998 


6 


15.9043 


0.001391 


1.87 X 10" 


-7 


1.0000000 



Table 5: Illustration of the speed of saturation in excitation number n of the total 
hadronic width T Hq by partial widths r„ from exclusive channels of mass M„. In 
this example, mQ = 15/3, 771^ = 10/3, and m = 0.56/5, for which M//^ = 16.3113/5. 



rriq 


N 


r/fo/roPE 


-1 




-1 


5.0 


18 


1.4 X 10" 


'8 


4.1 X 10- 


-4 


6.0 


16 


1.2 X 10" 


-7 


4.1 X 10- 


-4 


7.0 


13 


3.2 X 10" 


^7 


4.1 X 10- 


-4 


8.0 


11 


6.3 X 10- 


-7 


4.1 X 10- 


-4 


9.0 


8 


1.2 X 10- 


-6 


4.2 X 10- 


-4 


10.0 


6 


2.2 X 10- 


-6 


4.2 X 10- 


-4 


11.0 


4 


4.3 X 10" 


~6 


4.2 X 10- 


-4 


12.0 


3 


9.0 X 10- 


~6 


4.2 X 10- 


-4 


12.5 


2 


2.1 X 10- 


-5 


4.4 X 10- 


-4 


13.0 


1 


3.2 X 10- 


-5 


4.5 X 10- 


-4 


13.5 


1 


3.4 X 10" 


~5 


4.5 X 10- 


-4 


14.0 





5.9 X 10" 


-4 


1.0 X 10- 


-3 


14.5 





8.4 X 10" 


-4 


1.3 X 10- 


-3 



Table 6: Numbers relevant to local duality violation for niq = 15/5, m = 0.56/5, and 
niq variable. indicates the excitation number of the heaviest final-state meson 
kinematically allowed for the given initial meson mass M^q'- Mj^ < M^q < M^v+i. 
Thq is the total hadronic width Eq. (|37|) , while Tope and Tq are given by Eqs. (^) 
and (p9|), respectively. 
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AT 
iV 


i / J- OPE — J- 


F /F 1 


o.u 





D.y X iu 


y.4 X iU 


0.0 





o.O X iU 


y.4 X iU 


6.0 


4 


1.3 X 10-5 


9.5 X 10-^ 


6.5 


3 


1.8 X 10-5 


9.5 X 10-^ 


7.0 


3 


2.1 X 10-5 


9.5 X 10-^ 


7.5 


2 


5.5 X 10-5 


9.9 X 10-"^ 


8.0 


1 


7.6 X 10-5 


1.0 X 10-3 


8.5 


1 


7.9 X 10-5 


1.0 X 10-3 


9.0 





1.4 X 10-3 


2.4 X 10-3 


9.5 





1.9 X 10-3 


2.8 X 10-3 


iable /: 


bame 


as in Table except rriQ = 10/3. 


rrig 


AT 
iV 


^Hq/^OPE — ^ 


r^o/rg-i 


1.0 


3 


8.5 X 10-5 


3.7 X 10-3 


1.5 


3 


8.0 X 10-5 


3.7 X 10-3 


2.0 


2 


1.9 X 10-^ 


3.8 X 10-3 


2.5 


2 


2.5 X 10-^ 


3.8 X 10-3 


3.0 


1 


3.5 X 10-^ 


3.9 X 10-3 


3.5 


1 


3.2 X 10-^ 


3.9 X 10-3 


4.0 





6.6 X 10-3 


1.0 X 10-2 


4.5 





7.3 X 10-3 


1.1 X 10-2 


iabie 8 


: Same 


as in Table ^ except mg = 5/3. 




AT 
iV 


r/fg/roPE — 1 




i.O 





1.3 X 10-^ 


1.8 X 10-1 


z.U 





2.4 X 10-2 


4.2 X 10-2 


O.U 


i 


1.5 X 10-3 


1.1 X 10-2 


A n 
4.U 


Z 


2.7 X 10-^ 


5.7 X 10-3 


5.0 


4 


6.0 X 10-5 


3.6 X 10-3 


6.0 


5 


2.7 X 10-5 


2.6 X 10-3 


7.0 


7 


1.1 X 10-5 


1.9 X 10-3 


8.0 


8 


5.7 X 10-6 


1.5 X 10-3 


9.0 


10 


2.8 X 10-6 


1.2 X 10-3 


10.0 


13 


1.3 X 10-6 


9.3 X 10-^ 


11.0 


15 


4.5 X 10-^ 


7.7 X 10-^ 


12.0 


18 


4.8 X 10-^ 


6.5 X 10-^ 



Table 9: Same as in Table except mg = m = 0.56(3 fixed and mq variable. 
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N 


^Hq/^ OPE 


-1 




-1 


1.0 





2.2 X 10- 


-1 


3.1 X 10- 


-1 


2.0 


1 


1.2 X 10- 


-3 


2.9 X 10- 


-2 


3.0 


1 


1.1 X 10" 


-3 


1.5 X 10- 


-2 


4.0 


2 


1.9 X 10- 


-4 


8.4 X 10" 


-3 


5.0 


4 


5.2 X 10- 


-5 


5.5 X 10" 


-3 


6.0 


5 


2.4 X 10" 


-5 


3.8 X 10- 


-3 


7.0 


7 


1.0 X 10" 


-5 


2.8 X 10- 


-3 


8.0 


8 


5.1 X 10- 


-6 


2.2 X 10" 


-3 


9.0 


10 


2.5 X 10- 


-6 


1.7 X 10- 


-3 


10.0 


12 


1.1 X 10- 


-6 


1.4 X 10- 


-3 


11.0 


15 


4.2 X 10- 


-7 


1.1 X 10^ 


-3 


12.0 


17 


5.2 X 10- 


^8 


9.7 X 10- 


-4 



Table 10: Same as in Table except mg = m = 0.26/3. 



n 


M2n 


R{MiJ X 10^ 


i?o(M|Jxl04 





1.7792 


930.5 


805.6 


1 


4.5349 


15.18 


15.30 


2 


6.2796 


4.082 


4.099 


3 


7.6574 


1.838 


1.844 


4 


8.8310 


1.037 


1.040 


5 


9.8703 


0.6635 


0.6651 


6 


10.813 


0.4602 


0.4614 


7 


11.681 


0.3376 


0.3385 


8 


12.490 


0.2581 


0.2588 


9 


13.250 


0.2036 


0.2042 



Table 11: Saturation of the vacuum current correlator as depicted graphically in 
Fig. 16. The fiducial point in each interval M|„_^ < < M|^_,_]^, in which R (aver- 



aged hadronic 5-functions) and Rq (Born-term partonic expression) are compared, is 
chosen to be M|^. 
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oo 50 25 15 10 



Figure 1: Values of the ground-state energy Mhq — mq versus niQ for m = 0.56/3 
determined through direct numerical calculation (solid line); fit of the mass expan- 
sion Eq. (H) to 0(l/mQ), using the relations Eq. (P) and the approach described in 
Appendix y (dashed line). 
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oo 50 25 15 
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Figure 2: Values of the ground-state integral cq y/mg for m = 0.56/3 (solid line). 
Co is related to the decay constant via Eq. (p8|). The dashed line represents an 
approximation using the first expansion in Eq. ([1^), good to O^l/mg), and the 
asymptotic value F from a polynomial fit in l/mq to the points on the solid line. 
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1.35 I > > r 




Figure 3: Values of the ground-state matrix element mQ(l — a;) for m = 0.56/3 (solid 
line). The dashed line represents an approximation using the second expansion in 
Eq. (|14D, good to 0(l/m|). 




0.1 - 



I ' ' ' ' ' ^ ' 

oo 50 25 15 10 7 5 3 

Figure 4: Values of the ground-state matrix element m^{{x'^) —{x)"^) for m = 0.56/5 
(solid line). The dashed line represents an approximation using the third expansion 
in Eq. (|T|), good to ©(l/mg). 



38 



0.014 
0.012 



0.01 - 




Figure 5: Values of the ground-state matrix element {mg / Mhq) {I / x)—l for m = 0.56/9 
(solid line). The dashed line represents an approximation using the the final expansion 
in Eq. (0), good to 0(l/m^). 
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mQ 

Figure 6: Same as in Fig. 0, except for the first excited state. 
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o 



'o 1-01 




0.99 



0.98 
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Figure 7: Same as in Fig. |T|, except m = 0.26/3. 
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LU 
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m 

l-H 
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Figure 8: Duality deviation between exact hadronic width and F determined from 
the OPE [Eq. or (|T8|)], good to 0(l/m^). Here = 15/3, m = 0.56/?, and 
variable. 
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Figure 13: Same masses as in Fig. 11, using the Born-term partonic rate instead 
of Fqpe- Note that the deviation is much larger. 



1/rxdr/dE (mb=10) 
From (2/mQ) x 9^ (mQ=25) 



1.2 
1 

0.8 
0.6 
0.4 
0.2 


0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

2E/mB 

Figure 14: Exact differential width T'^ dT/dE, averaged as described in the text, 
compared to the continuous parton distribution computed via Eq. (|i3|) with rriQ = 25/5, 
at mh = lOP and m = 0.56/5. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



1-2E/Mb 

Figure 15: The partially-integrated (in E) differential distribution ^{1—2E/Mb) and 
the corresponding smeared exact result dF as defined in the text, for the same 

inputs as in Fig. 0. 




20 40 60 80 100 120 140 160 

Figure 16: Saturation of the leading term in the OPE of the vacuum polarization 
function Rq by exclusive channel (5-function contributions, smeared as described in 
the text (R). 
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